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Abstract 

This work analyzes the problem of community structure in real-world net- 
works based on the synchronization of nonidentical coupled chaotic Rossler 
oscillators each one characterized by a defined natural frequency, and coupled 
according to a predefined network topology. The interaction scheme contem- 
plates an uniformly increasing coupling force to simulate a society in which 
the association between the agents grows in time. To enhance the stability 
of the correlated states that could emerge from the synchronization process, 
we propose a parameterless mechanism that adapts the characteristic fre- 
quencies of coupled oscillators according to a dynamic connectivity matrix 
deduced from correlated data. We show that the characteristic frequency 
vector that results from the adaptation mechanism reveals the underlying 
community structure present in the network. 
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1. Introduction 

Community structure appears in many networked systems, including a 
variety of biological, social, technological, and information networks PJ. A 
network can be represented as a graph G composed by a large number of 
highly interconnected units, and a community is characterized by a large 
number of edges connecting vertices within individual groups, with only low 
concentrations of edges between these groups |2j]. The aim of community 
detection in graphs is to identify the modules and, possibly, their hierarchical 
organization, by only using the information encoded in the graph topology 
The problem has a long tradition and it has appeared in various forms in 
several disciplines 0, H,!!, 0]- The general notion of community structure in 
complex networks was first pointed out in the physics literature by Girvan 
and Newman [3]. During the last years many algorithms have been proposed 
ranging from traditional methods (Graph partitioning, spectral clustering), 
modularity-based methods and synchronization-based dynamics algorithms 
@, S, @| • For a recent review about the problem of detecting communities in 



graphs see [10 



Dynamics algorithms based on synchronization have been studied by 



several authors [11|, 1121, |13|, 1141 ]. and have shown that the dynamical syn- 
chronization process in complex networks is a suitable method to unravel 
their different topological scales. Most of these works are based on the Ku- 



ramoto model [15| or on its variants such as the opinion changing rate model 



16[, where the dynamics of each isolated oscillator is modeled in a one- 



dimensional state space. However, as actual systems are always described 
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by high-dimensional complex dynamics that involve many variables and pa- 
rameters, it is of practical interest to know about community structure in 
complex networks in these circumstances. This issue was approached in 18 
considering N identical Rossler oscillators coupled according to the dynamics 
of the opinion changing rate model, making use of the edge betweenness load 
matrix which implies an additional computational cost. On the other hand, 
these algorithms are usually dependent on initial conditions and frequencies 
distribution of coupled oscillators, and give rise to unstable solutions. 

In this work we model the dynamics of a given complex network using 
an ensemble of nonidentical coupled chaotic Rossler systems, each one char- 
acterized by a defined natural frequency. Such a model is chosen because 
Rossler system has been fairly well studied, in particular, in the context of 



complex networks [17|, [18|, |19|, |20|, [21| . Others oscillatory systems can also be 
considered. 

Our proposed approach tackles the problem of community structure in 
networked systems focusing on two aspects, the first one is to consider an 
ensemble of interacting nonidentical chaotic Rossler systems where the cou- 
pling force increases linearly with time simulating a society in which the 
interconnections between the agents grow in time. In fact, the increase of 
interactions between the members of a society helps to stabilize a plurality of 



different non-interacting clusters [161 ] . The values of coupling strength range 
from zero, which corresponds to an initial uncorrelated state, to a maxi- 
mum value given by the ratio A max /A2, representing an intermediate state 
of highly interconnected units forming local synchronized clusters. A2 and 
A m ax are the first non-zero and maximum eigenvalues of the Laplacian matrix 
respectively. The Laplacian matrix of a graph is given by L = D — A, where 
D is the diagonal degree matrix and A is the adjacency matrix of the graph 
G. The connection between the spectral information of the Laplacian matrix 
and the hierarchical process of emergence of communities has been reported 



111 



11 



The second aspect is related to the influence that neighbors have to change 
the natural frequency (or opinion) of a particular oscillator (or agent), and 
therefore the emergence of the modular structure found in the network where 
each module will correspond to oscillators with common average frequency. 
Initially, when the coupling is zero, each oscillator evolves according to its 
natural frequency giving rise to uncorrelated states. As soon as the coupling 
is engaged new correlated states emerge corresponding to oscillators that will 
evolve at the same average frequency. We show that a dynamical adaptation 



3 



of the characteristic frequencies of oscillators, according to the idea of confi- 



dence bound [22] , improves the stability of these emergent correlated states, 
and gives rise to a final frequency vector that brings into relief the network 
community structure. To this end, a parameterless frequency adaptation 
approach is proposed. 

The paper is organized as follows. In section 2, the dynamics of the 
Rossler oscillators network is introduced. In section 3 we will report the 
proposed parameterless and enhanced adaptive approach used to change the 
characteristic frequencies of coupled oscillators. In sections 4 and 5 we de- 
scribe both the social and computer-generated networks considered to test 
our approach. Finally, the results of numerical simulations are drawn in 
section 6. 



2. Network dynamics 

The dynamics of a graph G of N nonidentical coupled oscillators can be 
described as follow: 

3 

where i = 1,...,N. = F(^), & G 9ft n is the local dynamics of each 

autonomous oscillator. Each oscillator can be characterized by a defined 
natural frequency which is the long-time average of the phase velocity. The 
oscillators are coupled with the coupling coefficients over a predefined 
topology characterized by the connectivity matrix Ay-, indicates the 
coupling strength between the ith and jth oscillators. If oscillators % and 
j are connected, they have a coupling strength = K\ otherwise, the 
coupling strength is Kij = 0. 

In this work, we study an ensemble of many interacting Rossler systems 
coupled through the y component: 

±i(t) = -LOiVi - z h 

ij.it) = uJiXi + ay { + g J2f A iAVi ~ Vi)> ( 2 ) 
Zi(t) = b+ (xi- c)zi. 

All Rossler oscillators are nonidentical, and the parameter cjj = loq + Acj. 
selected randomly from a uniform distribution corresponds to the natural 
frequency of the individual oscillator [23]. Awj is the frequency mismatch 
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between neighboring chaotic oscillators. We set a = 0.2, b = 0.2, c = 5.7 and 
loq = 0.9 so as to ensure the individual oscillator generates chaotic dynamics 
with phase coherent attractor. k{ is the degree of the ith oscillator, K is the 
coupling strength, and Ay are the elements of the adjacency matrix of G 
(Ajj are if vertices i, j are not connected and 1 otherwise). The y-coupling 
would guarantee the stability of the synchronized state for a coupling force 
large enough 24 . 



3. The adaptive mechanism of characteristic frequencies 

For a coupling force K = 0, each oscillator evolves according to its own 
dynamics resulting in uncorrelated states space variables. As soon as the 
coupling strength increases, it may be that several oscillators synchronize 
and oscillate with a common average frequency, whereas their neighbors have 
their own, different, frequencies. To enhance the stability of these emergent 
regimes, we propose a parameter-free dynamical adaptation mechanism to 
change the characteristic frequencies of oscillators. 

We assume that each oscillator repeatedly changes its characteristic fre- 
quency by the median frequency of their oscillators of confidence, i.e, oscilla- 
tors whose dynamics are confined within the confidence manifold. To define 
this confidence manifold, we proceed as follows: throughout the synchronized 
process we monitor the temporal evolution of the y component of each os- 
cillator during a predefined time interval. Then, we compute the similarity 
matrix S between the resulting time series. This matrix will be based on nor- 
malized Pearson's coefficients, and calculated as follows: let p^ denotes the 
Pearson's coefficient between two arbitrary time series yi and yj correspond- 
ing to the temporal evolution of the y component of the oscillators i and j 
respectively. These Pearson's coefficients are transformed to similarities that 
belong to the interval [0, 1]: 

Sij = (Pij - min(pij))/(l - min(pij)) 

where min(pij) denotes the minimum p^ over the 5"s entries. Based on this 
similarity matrix, for each oscillator i, we can compute its average similarity 
to all other oscillators. Let m.j denote this average. 

The matrix S is converted to a dynamic connectivity matrix A* that will 
reflect the dynamic of real associations between the oscillators during the 
synchronization process. This matrix is constructed according to an adaptive 
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neighborhood process as follows: 

A* = { Ai i if Si i > mi or Si i > m i 
l i otherwise 

Finally, for each oscillator i, we define its neighbors of confidence, 
as the set of adjacent oscillators according to A*. So, the frequency of the 
oscillator i is replaced by the median frequency of oscillators forming the 
confidence neighbors. Using the median instead of the mean frequency tends 
to be more robust to outliers, and gives rise to a frequency vector with a 
clear clustering structure. 

This mechanism of frequencies adaptation is initiated once the coupled 
oscillators reach a certain degree of synchronization characterized by high 
values of average similarity m,. We first consider a phase in which the cou- 
pled oscillators are forced to synchronize through the increase of the cou- 
pling strength and according to the underlying network topology. When the 



coupling strength K reaches the appropriate value [20j for which complete 
synchronization is not established, but others synchronization features could 
emerge, we keep it fixed. Then, we began the mechanism of frequencies 
adaptation as described above. This mechanism takes place every few time 
units (or adaptation step) during a specific time interval. Both the simu- 
lation time and the adaptation step determine the length of the time series 
used to compute the similarity matrix. The number of times we adapt the 
oscillator frequencies is given by the ratio between the time interval and the 
adaptation step, and should be limited to avoid a homogeneous frequency 
vector. 

The adaptation procedure gives rise to a one-dimensional frequency vec- 
tor with some clustering structure able to detect the underlying modules 
present in a given network. Each adaptation step a new frequency vector is 
obtained, and then new community subdivisions will emerge. So to deter- 
mine the optimal number of modules or clusters we have adopted the criterion 



of maximum modularity [25|, [26| . The modularity is defined as the fraction 
of links within communities minus the expected fraction of such links in a 
random network. This measure provides a way to determine if a certain de- 
scription of the graph in terms of communities is more or less accurate. High 
values of modularity should indicate good partitions with many more inter- 
nal connections than expected at random. This criterion has been adopted 
for two reasons: is the most widely used in the literature, and the algorithms 
chosen to compare our methods adopt the same criterion. Although recently, 
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it has been shown that the modularity maximization could exhibit extreme 
degeneracies admitting an exponential number of distinct high-modularity 
solutions depending on both the size of the network and on the number of 



modules it contains [27J. In this work, the existence of possible degenerate 
solutions is limited since we have considered small real-world networks with 
few number of modules. 



4. Real- world networks 

In this paper, we have considered four real networks from the social sci- 
ence literature which have grown into standard tests for community structure 
detection algorithms. These networks are: 

- The karate club network |28| which shows the partition of friendships 
between the members of a karate club at US university in the 1970s. It 
consists of 34 members, and after the observation and construction of 
the network, the club in question splits in two as a result of an internal 
dispute. 

- The network of bottlenose dolphins [29| that was constructed from ob- 
servation of a community of 62 bottlenose dolphins, where ties between 
dolphin pairs were established by observation of statistically significant 
frequent association. Bottlenose dolphins communities have been de- 
scribed as a fission-fusion societies and therefore individuals (or agents) 
can make decisions to join or leave a group. The community structure 
of this network reveals two major communities. The first one is formed 
by 21 dolphins, while the second community which further splits into 
three sub-communities groups all other dolphins. 

- The jazz bands network [30] which consists of 198 bands performed 
between 1912 and 1942. Two bands are connected if they have a mu- 
sician in common. The bands network highlights the existence of two 
large communities, where the largest community also branches into two 
sub-communit ies . 

- The network of American football games between Division IA colleges 
during regular season Fall 2000. Nodes in the network represent teams 
and edges represent games between teams [2]. This network is char- 
acterized by a known community structure in which nodes are divided 
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into conferences containing around 8-12 teams each. Games are more 
frequent between members of the same conference than between mem- 
bers of different conferences. 



5. Tests on computer-generated networks 

Moreover, we have tested the performance of our method on computer- 
generated graphs with a known community structure. We have adopted 
the early class of the standard benchmark graphs of Girvan and Newman 



introduced in [31[. Theses graphs are generated with N = 128 nodes, split 
into four communities of 32 nodes each. Links between nodes belonging to 
same community are drawn with probability Pi n , while a pair which belong 
to different community are joined with probability p out . The value of p out is 
chosen so that the average number of inter-community edges per node, i.e, 
z out , can be controlled. Each node have on average Zi n edges to nodes in the 
same community and z out edges to nodes in other communities, maintaining 
a total average node degree k = Z{ n + z out = 16. 

We consider values of z out ranging from five, corresponding to clear com- 
munity structure (z out <C z in ), to 10 which describes a network with bad 
defined structure [z out 3> Zi n ). Since the "real" community structure is well 
known for these trial networks, we can validate our method computing the 
fraction of correctly identified nodes. 



6. Numerical results 

In this section, we integrate the system of Eq. [2] over the real-world 
networks described in the previous section. Each network is characterized 
by a connectivity matrix Aij, and the initial characteristic frequencies Ui 
which distinguish each oscillator. Values of Ui are randomly chosen from 
a uniform distribution with a frequency mismatch Aui = 0.2, so that the 
Rossler oscillators will evolve in a chaotic phase- coherent regime. 

Numerical results have been performed simulating the system over 75 
time units. We divided the simulation process into two phases. During the 
first stage (time < 20), every one time units, the coupling force is increased 
linearly from an initial value K = to a value not large enough to induce 
a regime close to complete synchronization. The values achieved are Kz = 
1.50, K D = 1.80, Kj = 2.13 and K F = 0.9793 for the Zachary, Dolphins, 
Jazz and Football networks respectively. In this stage, the coupled oscillators 
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are forced to synchronize mutually according to the network topology and 
the interaction strength. 




Time 



Figure 1: The average similarity of the oscillator networks versus time performed over 100 
runs of the algorithm. 

In the second phase (time >= 20), we keep K fixed and then we began the 
process of frequencies adaptation. As for the coupling force, every one time 
units, we adjust the oscillator frequencies according to the adaptation mech- 
anism and we register values of average similarity (m.j) of each oscillator to 
characterize the temporal coherence of the network, and values of modularity 
Q to quantify the goodness of a partition of the network into communities. 
To compute Q, the nodes in the network are clustered according to the re- 
sulting one- dimensional frequency vector. For all the networks considered in 
this work, this vector has shown a clear community structure corresponding 
to oscillators that evolve at the same average frequency. To determine the 
optimal number of communities, we have computed the modularity for a di- 
vision of a network into different number of communities ranging from two to 
eight communities. The values of modularity achieved by our approach have 
been compared with those obtained by the widely used betweenness-based 
algorithm of Girvan and Newman ji[ . 
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In Fig.l we report the global temporal coherence of the whole network 
[(rrii)}, where (.) stands for the average of all oscillators and [.] denotes aver- 
aging over 100 different system simulations. 




Figure 2: (a) Time evolution of the maximum modularity max{Q) achieved over a set 
of 100 different runs of the algorithm, (b) Maximum values of modularity corresponding 
to each run and reflecting the variability over the 100 different realizations. For Zachary 
and Dolphins the maximum modularity corresponds to a division into four communities. 
For the Jazz network the optimal division reveals three communities, while for Football 
network the optimal modularity is achieved for a partition into eight communities. In the 
inset of panel (b) temporal modularity of the Zachary network corresponding to one of 
the optimal runs. 

As it can be appreciated, during the first stage, the network coherence 
shows an increasing pattern due to the gradual increase of the coupling force. 
It is important to stress that, in this region of the coupling strength, agents 
in each final group tend to maintain a synchronized regime despite their 
different characteristic frequencies. This is due to the strong influence the 
community structure present in the network which affects notably the emer- 
gence of synchronized groups. After the adaptation mechanism takes place, 
we appreciate a slight increase of the global coherence of the network that 
we associate with the emergence of synchronized groups oscillating at the 
same average frequency. In this phase, as it can be seen in Fig 2 (a), we 
observe that the effect of the adaptation is immediate causing a sharp in- 
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Network 


Qn (N c ) 


Qmax (-^C ) 


size 


Zachary 

Dolphins 

Jazz 

Football 


0.4090 (4); 0.4010 (5) 
0.4580 (4); 0.5190 (5) 
0.4379 (4); 0.4452 (5) 

0.5470 (6); 0.5980 (8) 


0.4174 (4); 0.4143 (5) 
0.5220 (4);0.5184 (5) 
0.4437 (4); 0.4436 (5) 
0.4440 (3) 

0.5746 (6); 0.5901 (8) 


[13,11,5,5] 
[20,12,9,21] 

[61,75,62] 

[20,10,15,25,12,11,9,13] 



Table 1: Comparison of modularities for the network division found using our approach 
and betweenness-based algorithm of Girvan and Newman [Girvan et al. 2002]. Nc is 
the number of communities found. Qn is the modularity achieved by the algorithm of 
Girvan and Newman. Qmax is the maximum modularity corresponding to the optimal 
configuration. Also included is the size of communities corresponding to the optimal 
partition. The number of communities with maximum values of Q have been highlighted. 



crease in the values of the temporal modularity Q(t). In fact, Fig. 2(a) 
reports the maximum values of modularity Q ma x achieved over a set of 100 
different executions versus time. These values correspond to a division into 
four communities for the Zachary and Dolphins networks, while for the Jazz 
and Football networks the maximum values are achieved for a partition into 
three and eight communities respectively. All these subdivisions appear to 
show stable patterns with almost the same values of Q max in the case of the 
Dolphins and Jazz networks, and oscillating values with small deviations for 
the Zachary and football networks. On the other hand, according to the 
variability of the maximum values of modularity over the 100 runs of the 
algorithm (Fig. 2(b)), we can conclude that, in general, with few runs of the 
algorithm one can reach the maximum modularity values reported in Table 
1. Data reported in Fig. 2(b) corresponds to sorted values of maximum 
modularity over the set of 100 runs of the algorithm. Indeed, we can appre- 
ciate that for the Zachary network 83% of time, Q max > 0.363 which is the 
value achieved by the edge-betweenness algorithm for a partition into four 
communities. For the Dolphins network, Q m ax is greater to the value 0.458 
associated with a partition into four communities by the edge-betweenness 
algorithm with a percentage of 88%. The inset of Fig. 2(b) shows for the 
Zachary network an example of time evolution of modularity corresponding 
to one of the optimal runs. 

In Table 1 we report values of modularity corresponding to the optimal 
partitions achieved taking our approach compared to those obtained by the 
betweenness-based algorithm of Girvan and Newman jij. Depending on the 
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number of communities found using the two approaches, the difference in 
maximum modularity is about 5.5% for small number of communities and 
1% when the number of communities is somewhat higher. 

In the following we present, for each of the considered networks, the 
numerical outcomes illustrating the clustering of the characteristic frequency 
vector resulting from the adaptation process. 



6.1. The Zachary network 

In Fig. 3(a) we report values of Ui versus node index in the Zachary 
network for a partition into two communities in accord with the natural 



division originally observed in [28| . The maximum modularity value achieved 
by our approach for this partition is Q z = 0.3718. As it can be appreciated, 
Ui values are visibly split into two main groups. Nodes 3 and 10 are assigned 
to the erroneous community and display values of Ui slightly lower than the 
average frequency (a;) 1 = 0.9486 associated with the community to which 
they belong. Fig. 3(b) shows values of Ui sorted in ascending order (circle 
marker) contrasted with the average frequency of each community (dashed 
line). 

On the other hand, the maximum modularity achieved by our approach 
is reached for a division of the Zachary network into four communities with 
a value Q^ ax = 0.4174 which is 14% higher than the modularity obtained 
by the betweenness-based algorithm for the same partition. The resulting 
partition represents the division of the two initial communities into two sub- 
communities each one. The first community consists of 13 nodes oscillating 
at the average frequency (u) 1 = 0.8870, while the second one with 11 nodes 
is characterized by the average frequency (u) 2 = 0.8949. The last two com- 
munities both with five nodes each are identified by the average frequencies 
(u) 3 = 0.8958 and (w) 4 = 0.8976 respectively. 



6.2. The Dolphins network 

The natural division of the Dolphins network originally observed in 29 
shows two major communities where one of them is further divided into 
three sub-communities. This partition has been detected using our adapta- 
tion mechanism leading to a frequency vector whose values versus node 
index are depicted in Fig. 4(a-b). As it can seen, we can associate the first 
community to oscillators whose average frequency is about (u) 1 = 0.8768. 
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Figure 3: The Zachary network, (a) The scatter plot illustrating the clustering of the 
characteristic frequencies of oscillators, (b) The sorted values of frequency of oscillators 
(circle marker) according to a division into two clusters achieved by our approach. In 
dashed line is represented the average frequency of each community. 



The number of nodes in this community is 21. All other oscillators are part of 
the second community which it could be divided into three sub-communities 
oscillating at the average frequencies (u) 2 = 0.9497, (u) 3 = 0.9617 and 
{u) A = 0.9433 respectively The distribution of nodes in these three sub- 
communities is as follow: 20 oscillators move at the average frequency (w) 2 , 
12 at (u) 3 and 9 at (u)) A . The modularity corresponding to this division has 
a value Q D = 0.5220 which is about 13% higher than the value achieved by 
the betweenness-based algorithm. 



6.3. The Jazz network 

In Fig. 5(a) we report values of Ui versus node index in the Jazz network 
for a partition into three communities in accord with the natural division 
originally observed in 30j ■ The maximum modularity value achieved by our 
approach for this partition is Q J = 0.444. As it can be appreciated, the 
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Figure 4: The Dolphins network, (a) The scatter plot illustrating the clustering of the 
characteristic frequencies of oscillators, (b) The sorted values of frequency of oscillators 
(circle marker) according to a division into four clusters achieved by our approach. In 
dashed line is represented the average frequency of each community. 



three communities are clearly identified by the split of the frequency vector. 
The first one characterized by the average frequency (u) 1 = 0.8765 is made 
up of 61 nodes. The second community with 75 nodes can be recognized by 
the characteristic average frequency (a;) 2 = 0.9232. In this case, we observe 
two nodes (nodes 17 and 74) that can be considered as outliers since they 
oscillate at a frequency much higher than the frequency characterizing this 
community. Finally, the third community is formed by 62 nodes oscillating 
at the average frequency (u) 3 = 0.9705. 



6.4- The Football network 

Finally, as for the previous networks, the Fig. 6(a-b) report values of Ui 
versus node index in the Football network for a partition into eight commu- 
nities. One can clearly observe how the adaptive mechanism clusterizes the 
characteristic frequency vector resulting in oscillators that evolve at the same 
average frequency. The maximum modularity value achieved by our approach 
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Figure 5: The Jazz network, (a) The scatter plot illustrating the clustering of the charac- 
teristic frequencies of oscillators, (b) The sorted values of frequency of oscillators (circle 
marker) according to a division into three clusters achieved by our approach. In dashed 
line is represented the average frequency of each community. 

for this partition is Q F = 0.5901 which agree with the value achieved by the 
betweenness-based algorithm. This partition in eight communities is charac- 
terized by the following average frequencies: (u) 1 = 0.8428, (u) 2 = 0.8615, 
(u) 3 = 0.8767, (w) 4 = 0.8949, (u) 5 = 0.9244, (u) e = 0.9454, (u) 7 = 0.9490 
and (u) s = 0.9508. And corresponds respectively to clusters with 12, 20, 25, 
13, 9, 15, 10 and 11 nodes each. Further clustering of some identified com- 
munities, gives away that the oscillators with the average frequencies (a;) 2 
(square marker) and (u>) 3 (circle marker) could be further divided into two 
sub- communities each one giving rise to a partition into ten communities. 



6.5. The ad hoc trial networks 

To conclude, we report experiment results carried out with ad hoc trial 
networks as described in section 5. Figure 7 shows values of the fraction 
of correctly identified nodes, averaged over twenty-five different realizations 
of the computer generated networks, as a function of the average number of 
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Figure 6: The football network, (a) The scatter plot illustrating the clustering of the 
characteristic frequencies of oscillators, (b) The sorted values of frequency of oscillators 
(circle marker) according to a division into eight clusters achieved by our approach. In 
dashed line is represented the average frequency of each community. 



inter-community edges z out . We have also registered values of the modularity 
Q achieved by our approach (see inset of Fig. 7). These results have been 
compared with those achieved by the synchronization-based dynamical clus- 
tering algorithm (OCR-HK) of [lij. The improvement in the performance 
of our approach is clearly appreciated as soon as z „ ut exceeds the value of 
six inter- modular edges per node (see Fig.l in Ref. [12|). The sensitivity of 
the OCR-HK algorithm starts to decay above z out > 6, while for our method 
this decrease starts for z out > 7. On the other hand, when the networks have 
a clear community structure (z ou t < 6) the performance of our method is 
complete achieving 100% of prediction correctly. For this value, the OCR- 
HK identify about 90% of the nodes correctly, while our method is able to 
correctly classify 94% of the nodes. For values of Z out > 7, both algorithms 
perform the same fraction of correctly classified nodes. It should be stressed 
that for values of Z out = 8 only few algorithms are still able to identify over 
80% of the nodes correctly (32[, but with a hight computational cost. 
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Figure 7: Fraction of the correctly identified nodes versus the number of inter-community 
edges per node Z out averaged over 25 graph realizations. In the inset we report values of 
modularity Q achieved by our approach (circle marker) and those corresponding to the 
real partitions (square markers) versus Z out . 

7. Conclusions 

This paper considers the problem of community structure in complex real- 
world networks based on the synchronization of nonidentical chaotic oscilla- 
tors. We have discussed a parameterless approach that adapts dynamically 
the characteristic frequencies of coupled oscillators to enhance the stability 
of the emergent synchronized groups according to a dynamic connectivity 
matrix. This adaptation mechanism gives rise to a frequency vector able 
to reveal the community structure presents in a given network. The pro- 
posed method has been tested both in real-world and ad hoc trial networks. 
The goodness of a given division of the network into communities has been 
quantified using the criterion of maximum modularity. The obtained values 
of modularity have been compared with those achieved by the betweenness- 
based algorithm of Girvan and Newman. The results are promising and fall 
within the scope of the results published in the literature for small networks 
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with few number of communities. It should be stressed that the approach 
proposed in this work was inspired by the synchronization-based dynamical 
clustering algorithm (OCR-HK) of 12J which scales as 0(N 2 ) tacking into 
account that it requires to compute the edge-betweenness of the network. So, 
in the worst case, our algorithm will scale as 0(N 2 ) too, which is a reasonable 
computational cost to address larger networks. 
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